On the role of step-flow advection during electromigration-induced step bunching 
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We propose a one-dimensional model based on the Burton-Cabrera-Frank equations to describe 
the electromigration-induced step bunching instability on vicinal surfaces. The step drift resulting 
from atomic evaporation and/or deposition is explicitly included in our model. A linear stability 
analysis reveals several stability inversions as the evaporation rate varies, while a deposition flux is 
shown to have a stabilizing effect. 
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Due to its importance for both fundamental science 
and technological applications in microelectronics, re- 
search in the field of crystal growth on semiconductor 
vicinal surfaces has become increasingly active 0, Q, H, 
S 0, El- Vicinal surfaces are slightly misoriented with 
respect to closed-packed crystalline planes of low indices, 
such as the well-studied Si(lll) plane. The atomic steps 
resulting from the miscut are generally mobile (step flow) 
and subject to a number of instabilities as step bunch- 
ing or meandering [H, 0, H, 0, [l(| Eli- At long times, 
nonlinear dynamics lead to the formation of more or 
less well-organized surface patterns. The idea to control 
these patterns is appealing, a major issue being the self- 
organization of nano-structures on semiconductor sur- 
faces. A few ways to reach this goal are currently under 
investigation, such as imposing a net atomic flux (evap- 
oration/deposition), using elastic stresses or applying a 
constant electrical field inducing an adatom drift (surface 
electromigration). Combinations of these methods allow 
extra flexibility [g, [f| • Surface electromigration was first 
observed by Latyshev and co-authors [13j |. and an early 
theory proposed by Stoyanov [ljj]. The stability of a 
Si(lll) vicinal surface may change according to the tem- 
perature, the sign of the net atomic flux, and the direc- 
tion of the electrical current. Three temperature regimes 
have been identified experimentally 15 , 3, 17 , 18| . Both 
in the high and low temperature regimes, the step bunch- 
ing instability appears when the electrical current is ap- 
plied in the step-down direction. Alternatively, in the 
intermediate regime, it may be necessary (according to 
the experimental conditions) to apply a step-up current 
to trigger step bunching. This stability inversion with 
respect to the direction of the electrical current is still 
the subject of active research and different mechanisms 
have been proposed such as sign variations of the effective 
charge number (Z*) with temperature, step transparency 
Hi , 2oj ] , and two-region terrace diffusion [2l| . While sign 



variations of Z* with temperature have been ruled out 
by a recent experiment of mass transport across a trench 
|22j |. the two other possibilities remain. However, due to 
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FIG. 1: Schematic side view of a vicinal surface: notations 
and definitions. The electrical field is imposed in the step- 
down direction here. 



the large number of experimental parameters and ultra- 
high vacuum conditions, direct in-situ experimental mea- 
surements of adatom diffusion at the step edges are still 
difficult to achieve. In this Letter, we propose a one- 
dimensional model based on the Burton-Cabrera-Frank 
equations [23| to describe the step flow instabilities aris- 
ing when electromigration is combined with a net atomic 
flux. We first remark that, under typical experimental 
conditions, the advection effects due to the mean step 
velocity are comparable in magnitude with the drift due 
to the electromigration current. Both effects are thus in- 
cluded in our model and a linear stability analysis shows 
that the interplay between them does provoke stability 
inversions as evaporation is increased. 

During evaporation or growth, the steps of Si(lll) vic- 
inal surfaces undergo a global drift with a constant mean 
velocity Vo- In the frame moving at this velocity, the 
adatom concentration C n on the n-th terrace obeys the 
following quasi-static equation, 
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(1) 



where D s is the diffusion coefficient, — kgT/ (Z*eE) 
the electrical length, r the desorption time, and F the 
atomic deposition flux. As shown in Fig. (fTJ), the electri- 
cal field reads E = EX., with X the unit vector pointing 
in the step-down direction. At terrace ends X — X n and 
X = X n+ i, Eq. (fT]) is subject to the following boundary 
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conditions, 
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D s /e E )c n = v[c n -ci% 

D s /t E )C n = -u[C n -C: q 



(2) 
(3) 



which ensure mass conservation. Here v is the step ki- 
netic coefficient, the adatom equilibrium concentra- 
tion at step n, and we neglect step transparency. We can 
remark from the structure of Eqs. |T][3]) that the advec- 
tion and the electromigration terms play a similar role. 
The associated velocities Vo and D s /£e may be both pos- 
itive or negative, according to the sign of the imposed 
net atomic flux and electrical current. The advection 
terms are usually neglected in theoretical studies, which 
is equivalent to setting Vo to zero in Eqs. ([US])- How- 
ever, on the basis of the present experimental knowledge, 
it seems unrealistic to neglect advection as compared 
to electromigration. On one hand, the net atomic flux 
varies in a rather wide range in practice (three decades) 
while the mean terrace width extends over one decade 
at least. As a consequence, the corresponding values of 
the mean step velocity Vq typically range from 10~ 10 to 
1CP 6 m.s -1 . On the other hand, due to the experimental 
uncertainties on the values of the effective charge number 
Z* and, to a less extent, to the temperature influence on 
the adatom diffusion coefficient D s , estimates of the drift 
velocity D s /£e are found to lie in the very same range. 
Excepted in equilibrium conditions, the two effects are 
thus similar both in nature and in magnitude, so that we 
will keep the Vo terms in our model. At a given time t, the 
velocity of step n is V n = Q s v[C n + C„_i — 2C%*] — Vq, 
where the concentrations are evaluated at X — X n (t), 
and f2 s represents the unit atomic surface. The terrace 
lengths L n (t) = X n+ i(t) — X n {t) vary slowly in time, 
with a constant mean value Lq. We introduce standard 
step-step repulsive interactions in the adatom equilib- 
rium concentration, C^/Cq = 1 + A(l/L n - 1/Lf l _ 1 ), 
where ksTAVl^ 1 is the step-step interaction coefficient 
[24] ]. Setting the unit length to the initial terrace width 
Lq and the unit time to to — Lq/ D s , we define the dimcn- 
sionless variables x = X/Lq, vq = VoIq/Lq, c n = C' n /Co, 
and c^ 9 = C^/Cq. For these variables, the previous 
equations become 

(d x + Vq - 1])d x C n - s 2 c n + f = 0, X n < X < X n +1, (4) 

x = x ni (5) 
x = x n+1 , (6) 
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with cf g = 1 + a(l/ln - 1/ln-i)- The system dynam- 
ics is thus controlled by six independent nondimcnsional 
parameters: 77 = Lq/£e is proportional to the electri- 
cal field, s 2 = Lq/(D s t) involves the rate of desorp- 
tion and / = FL^/(D s Cq) the atomic deposition flux, 
p = vLq/ ' D s compares attachment to diffusion, (f> = fl s Co 
gives the proportion of occupied sites, and a = A/Lq 




FIG. 2: Growth rate a as a function of the wave number fc, 
as obtained from Eqs. (|12l I13[) . Different evaporation rates 
s are considered. Nondimensional parameters are <j> — 0.2, 
/ = 0, 77 = -10~ 4 , p = 20.0, a = 5 x 10" 6 . 



measures step-step repulsion. These six parameters are 
not completely independent in practice. As a conse- 
quence, when temperature is varied in a given experi- 
mental setup, the system is driven along a rather com- 
plex trajectory in the parameter space. This trajectory 
is very likely to cross a succession of stable and unstable 
regions, as often observed in practice. 

To provide quantitative support to this idea, we now 
perform the linear stability analysis of a uniform train of 
steps traveling at a constant velocity in the laboratory 
frame. Let us remark that we keep the notation t for 
the dimcnsionless time hereafter. A general solution of 
Eq. (TjQ) giving the concentration profile on terrace n is 
Cn = f/s 2 + a n e r ^+/3 n e r2X , where n, 2 = -(v -rj)/2± 
1/2[(vq — ri) 2 + As 2 ] 1 / 2 . To avoid lengthy expressions, we 
first set / = and we discuss the case / > later on. 
For the unperturbed system, we find that vq obeys the 
following transcendental equation 



vo _ 2s 2 (u 1 - u 2 ) - puiU 2 A 
4>p d\ - d 2 



(8) 



where A = n — r%, Ui = e'' 1 — 1, U2 = e T2 — 1, 
di = (n + p)(r 2 - p)e Tl , and d 2 = (n - p)(r 2 + p)e V2 . 
The linear stability of the system is tested by adding a 
small harmonic perturbation of wave number k, on the 
step positions, x n = n + ee tkn+ ^ t . In order to find the 
dispersion relation £ = a + iuj, we expand the step veloc- 
ities up to the first order in the perturbation amplitude, 
e«l. We obtain the following exact expression for the 
growth rate: 



a(fc) = 20psin 2 (-) 



2 fk\ A 1 s 2 + a(A 2 + A 3 cos k) 



2 (di - d 2 ) 2 



(9) 



where 

Ai 
Ap 

A 2 = 6s 2 (d 2 -d 1 )[2s 2 {u 2 -u 1 )-pA(u 1 +u 2 + 2)], 

A 3 = 6s 2 pA(d 2 -d 1 )(l + e ri+r2 ). (10) 



[rid 2 - r 2 d x + 2pA( ri + r 2 )]e ri+r2 + r 2 d 2 - ndi, 
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FIG. 3: Stability diagram in the (s,rj) plane. Solid (dotted) 
curves represent rj c (rj p ). Nondimensional parameters are <j> = 
0.2, / = 0, p = 20.0, a = 10~ 4 . For a step-up electrical 
field, the horizontal dashed line represents a trajectory as the 
evaporation rate increases. 



In the parameter space, this allows us to compute the 
critical values of the parameters denning the boundaries 
between stable and unstable regions. One thus has to 
solve the system of two equations, Eqs. (JSJOi m the two 
unknowns, rj and vq. Remarking that the condition \vq — 
77I s is always verified for the practical values of the 
electrical field and net atomic flux, we obtain analytical 
expressions for the velocity, 
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p(l + e s ) +s(e s - 1)' 



and the growth rate, 
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with e? 3 = (s + p) 2 e — (s — p) , and 



-^- = ( s + p) 2 ( s -l)e 2 
ps e 

-,2s 



Aps e s + (s - pf(s + 1), 



B 2 = -6s 2 d 3 [(s + p)e 2s + p-s],B 3 = I2s 2 d 3 p e s .(13) 

Depending on the physical parameters, the growth rate 
a(k) assumes negative or positive values. In the latter 
case, the one-dimensional train of steps is linearly unsta- 
ble. The variations of the growth rate with the perturba- 
tion wave number k are displayed in Fig. ([2]) for differ- 
ent values of the evaporation parameter s. The electrical 
field is imposed in the step-up direction here (rj < 0). 
For s = 0.001 all the modes are stable, while the large 
wavelengths are unstable for s — 0.04. At higher evap- 
oration rates, all the wave numbers are unstable. Note 
that for s — 0.09, the most unstable mode is obtained for 
k = 7r (step-pairing instability). Using Eqs. (|11|12|13|) . 
we obtain the analytical expression for the boundary lines 
r) c {s) separating the stable and the step bunching regions 
in the (s, rj) plane: 



Vc 



£> 2 + B 3 
Bi 



(14) 



FIG. 4: Critical value of the nondimensional evaporation rate 
s* as a function of the attachment /diffusion parameter p for 



/ = (log- log scales). 



In addition, one shows that step-pairing instability be- 
comes the most unstable mode for r\ p — vq + a(B 2 — 
3B 3 )/Bi. The stability diagram representing the bound- 
aries between the stable, and unstable regions is shown 
in Fig. ((3|). For typical values of the step- up elec- 
trical field (500V/m), 77 assumes small negative values 
(lower than 10~ 2 in magnitude) which lie slightly be- 
low the horizontal axis in the (s,rj) plane [281 ] - Since 
the evaporation rate s becomes larger as temperature 
is increased, the system is stable both at low and high 
temperatures and it becomes unstable in the intermedi- 
ate temperature range. Indeed, this succession of stabil- 
ity reversals is observed experimentally for step-up cur- 
rents 15, 3, 17, [HI]. In the unstable regions, both for 
step-up and step-down currents, step bunching is super- 
seded by the step-pairing instability for large electrical 
fields. Moreover, the bunching regions progressively re- 
duce in size as the step-step interaction coefficient a is 
decreased. As shown by Eqs. (fl~3l[T4|) . the location of the 
vertical asymptote s — s* solely varies with the attach- 
ment/diffusion parameter p. These variations are repre- 
sented in Fig. ((¥]). A maximum is obtained for p ~ 1, 
when the attachment length d = D s jv is comparable to 
the terrace width Lq. For attachment-limited dynamics, 
p < 1, we find that s* ~ (12/9) 1 / 4 , while s* ~ (12/p) 1 / 2 
for diffusion- limited dynamics, p > 1. As discussed 
above, \rj\ < 10~ 2 in practice, so that the higher critical 
value of the evaporation rate is very close to s = s* . Since 
the desorption time reads r = Lq/(s 2 D s ), we respectively 
obtain t* = {Ll/UvD,) 1 / 2 , and t* = vLl/\2D 2 in 
the attachment-limited and the diffusion- limited regimes. 
This implies that s* strongly depends on the miscut an- 
gle. Let us discuss now the case of a nonzero deposition 
flux, / > 0. The drift velocity reads 



-20p 



Qs 2 -/)(e s -l) 
sp(l + e s ) + s 2 (e s - 1)' 



(15) 



It vanishes when the deposition and evaporation fluxes 
compensate, f = fo = s 2 . We first consider the case of 
a step- up current. The step bunching instability occurs 
only when evaporation dominates deposition, f < fo, as 
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FIG. 5: Stability diagrams in the (s, f—fo) plane. Horizontal 
axes f = fo correspond to zero net flux. On the left (right) 
side of the vertical axes, the electrical current is step-up (step- 
down). Solid lines represent the boundaries between stable 
and unstable regions. The regions below the dotted lines are 
unphysical (/ < 0). Nondimensional parameters are <j) = 0.2, 
p = 1.0, a = 10~ 4 , and we set a) r] = ±5.0 x 10~ 4 , b) r) = 
±5.0 x 10" 3 . 



shown on the left side of Fig. ([5]). On the other hand, for 
a step-down current, the instability may occur both for 
f < fa and / > /o, depending on the values of the phys- 
ical parameters. From Fig. ([5j) , it is also clear that in- 
creasing the deposition flux has a stabilizing effect, what- 
ever the direction of the electrical field. Since a ~ L^ a , 
after Eq. ([Li]) the stability threshold r\ c ~ d for large L 
values. Thus, measurements of the step velocity Vq could 
provide direct estimates of the effective charge number, 
Z* = (V k B T)/(D s eE). 

In conclusion, we have shown that advection can not 
in general be neglected as compared to electromigration. 
We have identified a set of six nondimensional param- 
eters governing the system evolution. For step-up elec- 
trical currents, two stability inversions have been found 
as evaporation increases. Numerical simulations of our 
model equations are currently performed to investigate 
nonlinear regimes, such as coarsening dynamics, oscillat- 
ing modes and spatiotemporal chaos 25, 2(| 27j ]. 
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